SLAC-PUB-14382 
ANL-HEP-PR-11-13 
NUHEP-TH/11-014 
BONN-TH-2011-05 



Supersymmetry Without Prejudice at the 7 TeV LHC 

John A. Conley 

Physikalisches Institut, Universitdt Bonn, Nuflallee 12, 53113 Bonn, Germany 

James S. Gainer 

High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439 USA and 
Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208 USA 

JoAnne L. Hewett, My Phuong Le, Thomas G. Rizzo 

SLAC National Accelerator Laboratory, 2575 Sand Hill Rd., Menlo Park, CA 94025, USA 

Abstract 

We investigate the model independent nature of the Supersymmetry search strate- 
gies at the 7 TeV LHC. To this end, we study the missing-transverse-energy-based 
searches developed by the ATLAS Collaboration that were essentially designed for 
mSUGRA. We simulate the signals for ~ 71k models in the 19-dimensional parameter 
space of the pMSSM. These models have been found to satisfy existing experimental 
and theoretical constraints and provide insight into general features of the MSSM with- 
out reference to a particular SUSY breaking scenario or any other assumptions at the 
GUT scale. Using backgrounds generated by ATLAS, we find that imprecise knowledge 
of these estimated backgrounds is a limiting factor in the potential discovery of these 
models and that some channels become systematics-limited at larger luminosities. As 
this systematic error is varied between 20-100%, roughly half to 90% of this model 
sample is observable with significance S > 5 for 1 fb _1 of integrated luminosity. We 
then examine the model characteristics for the cases which cannot be discovered and 
find several contributing factors. We find that a blanket statement that squarks and 
gluinos are excluded with masses below a specific value cannot be made. We next ex- 
plore possible modifications to the kinematic cuts in these analyses that may improve 
the pMSSM model coverage. Lastly, we examine the implications of a null search at 
the 7 TeV LHC in terms of the degree of fine-tuning that would be present in this 
model set and for sparticle production at the 500 GeV and 1 TeV Linear Collider. 



1 Introduction and Background 



The LHC has had an initial run at 7 TeV with both the ATLAS and CMS experiments 
collecting ~ 35 — 45 pb _1 of useful data. Even with this low integrated luminosity these 
experiments have been able to extend searches far beyond the reach of the Tevatron for many 
new physics scenarios with, so far, null results [HE]- This clearly demonstrates the power 
of increasing the center of mass energy in the search for new physics at hadron colliders. 
Starting soon, the LHC is to begin a longer run at 7 TeV and is expected to collect of order 
1 — 7 fb _1 of data over the next 2 years. Such a data set will allow for a first exploration of 
the TeV mass scale, and if new strongly interacting particles exist in this kinematic regime, 
they should be observed. 

A well-motivated, and perhaps most popular, possibility for new physics that may 
be discovered during this coming LHC run is Supersymmetry (SUSY) j3j. Both ATLAS [I] 
and CMS [5] have designed detailed searches for many of the SUSY partners of the Standard 
Model (SM) particles; these are mostly (but not exclusively) based on the assumption of 
mSUGRA/CMSSM-like soft breaking within the Minimal Supersymmetric Standard Model 
(MSSM) framework. This assumption greatly simplifies the exploration of the vast Super- 
symmetric parameter space. While these searches are designed to well cover the parameter 
space of these SUSY-breaking scenarios it is important to ascertain their discovery potential 
in a more general MSSM context. This is particularly worrisome in light of results from the 
Tevatron, where it has been realized IHJE] that relatively light gluinos may have escaped un- 
detected. The question then arises whether these specific scenarios and associated searches 
adequately describe the true breadth of the MSSM and its possible collider signatures, and 
whether the LHC searches as presently designed could fail to observe sparticle production. 
This has prompted several studies of more model independent search strategies as well as 
the development of simplified models [8]. 

Recently, we have addressed [9] this question by investigating the capability of the 14 
TeV LHC to explore a more general MSSM model parameter space, i.e., that of the pMSSM 
(phenomenological MSSM) [7J, to be described below, from the point of view of the ATLAS 
detector. In particular, we examined the performance of the planned ATLAS SUSY searches 
in exploring this more general MSSM scenario. This analysis provides insight into general 
features of the MSSM without reference to a particular SUSY breaking scenario or any other 
assumptions at the GUT scale. We found that the ATLAS mSUGRA-inspired searches, 
based on missing £r, did surprisingly well at covering the kinematically accessible portions 
of this model space; we also found that some interesting exceptions can arise in these more 
general models. Given the lower-energy run of the LHC over the next 2 years it behooves us 
to determine how well the corresponding mSUGRA motivated searches designed by ATLAS 
would perform under these conditions, since this is the situation presently before us. This is 
the goal of the present paper. We note that recently there have been several mSUGRA-based 
studies evaluating the capability of the 7 TeV LHC run to probe that parameter space [10] 
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As is well known, soft SUSY breaking within the MSSM in all generality leads to 
a theory with over ^ 100 a priori free parameters which prohibits a detailed study of this 
theory. A number of theoretically possible scenarios exist which describe the breaking of Su- 
persymmetry; maybe even multiple mechanisms are simultaneously responsible. Practically 
speaking, there are two ways to approach reducing this large number of a priori unknown 
parameters to something more manageable. One approach is to consider only specific, well- 
motivated SUSY breaking scenarios, such as mSUGRA or others. This leads to a drastic 
reduction in the number of free parameters (to only ^ 3 — 5) so that detailed analyses of the 
resulting parameter space can be easily achieved. A problem with performing such studies 
is that they may result in a bias as to the nature of SUSY signals when searching for collider 
or other SUSY signatures. An alternative approach is to be less prejudicial and to instead 
follow a bottom-up analysis which we have employed in a number of recent works P, [HJ and 
will make use of here. By imposing a set of theoretically and experimentally well-motivated 
constraints on the general MSSM (to be described below), without making any reference 
to the specific mechanism of SUSY breaking, we arrive at a theory with 19 TeV-scale pa- 
rameters. This is known as the pMSSM, which is significantly more manageable than the 
full Supersymmetric parameter space, and yet allows for more breadth than is present in , 
e.g., mSUGRAj^] These parameters will then completely define and describe all aspects of 
TeV-scale SUSY phenomenology. Such an approach has the advantage of being more general 
than any given (or given set of) specific SUSY breaking scenario(s) and allows one to be in 
some sense agnostic about the SUSY mass spectrum. 

To this end, we examine the £j? lss -based SUSY searches developed by the ATLAS 
collaboration for the 7 TeV LHC [12]. We simulate the pMSSM signal for roughly 71k 
pMSSM models (hereafter 'model' refers to a point in the 19-dimensional pMSSM parameter 
space) that we generated in our previous work [7j. We employ SM backgrounds provided by 
the ATLAS Collaboration. In Section 2, we describe our SUSY model generation and LHC 
analysis procedure. It is important to note that we strictly adhere to the analyses as designed 
by ATLAS. While numerous, and perhaps improved, SUSY collider search techniques have 
been discussed in the literature [13] , it is not our present purpose here to discuss or employ 
them. Section 3 contains our main results. We find the systematic error in determining 
the SM background to SUSY production is a limiting factor in the potential discovery of 
these models; in fact, some channels become systematics limited at larger luminosities. In 
this section we determine the fraction of our pMSSM model set that is discoverable at the 
7 TeV LHC. We then examine the model characteristics in some detail that render some of 
the models undetectable. We find that the observability of models depends on the precise 
details of the sparticle spectrum and that a blanket statement of constraints on the mass 
of, say the gluino or squarks, cannot be made. In this Section, we also explore potential 
modifications to the kinematic cuts in these analyses that may improve model observability. 
In Section 4, we examine the implications of a null SUSY search during this run with respect 

Even in such a case a full exploration of this large parameter space is at best difficult if not impossible 
with present computing power. 
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to the degree of fine-tuning present in these models, as well as the implications for sparticle 
production at a high energy e + e~ Linear Collider. A summary and our conclusions can be 
found in Section 5. 

2 Analysis Procedure for Inclusive SUSY Production 
at the LHC 

The purpose of this work is to explore how well mSUGRA-inspired inclusive SUSY searches 
(in particular the set proposed by the ATLAS collaboration [12J) apply to the larger and much 
more general pMSSM parameter space for the 7 TeV run of the LHC. This is similar in spirit 
to to Ref. [9], which explored this question for the more powerful 14 TeV LHC. The pMSSM 
model sample that we study was generated in Ref. [7j; here, we briefly review the procedure 
employed to generate this sample. We then describe our procedure for generating the signal 
events, comparing to background, and determining the statistical criteria for discovery. We 
will show that we faithfully reproduce the ATLAS results in each analysis channel for their 
benchmark SUSY model. 

As stated above, we study the 19-dimensional parameter space of the pMSSM. This 
set of parameters was arrived at [14] by imposing the following set of requirements onto the 
general R-Parity conserving MSSM: (i) the soft parameters are taken to be real so that there 
are no new CP-violating sources beyond those in the usual CKM matrix; (ii) Minimal Flavor 
Violation(MFV) [15] is assumed to be valid at the TeV scale; (Hi) the first two generations 
of sfermions having the same quantum numbers are taken to be degenerate and to have 
negligible Yukawa couplings and (iv) the lightest neutralino is assumed to be the Lightest 
Supersymmetric Particle(LSP) and is a stable thermal WIMP. Most of these assumptions 
are applied in order to avoid issues associated with flavor physics constraints. With these 
conditions, the remaining 19 free soft-breaking parameters are given by the three gaugino 
masses, M^=i_ 3 , the ten sfermion masses mj, the three A-terms for the third generation 
fermions (A^t), an d the usual Higgs sector parameters /i, Ma and tan/3. 

To generate the specific pMSSM parameter points that we study below (hereafter 
referred to as our set of models), we performed numerical scans over the space formed by 
these 19 parameters. This required both a selection of the parameter range intervals as well 
as an assumption about the nature of the scan prior for how points are chosen within these 
intervals. These issues are both described in detail in our previous works [9l [H] and the 
interested reader should refer to them directly. Here, we simply note that two scans were 
performed: one employing a flat prior beginning with 10 7 points and one with a log prior 
employing 2 x 10 6 points. The main distinctions between these two scans directly relevant 
to our analysis here are that (a) all SUSY mass parameters were restricted to be < 1 TeV 
for the flat prior case, while for the log case the upper limit on mass parameters was raised 
to 3 TeV, and (b) the choice of the log prior generally leads to a more compressed sparticle 
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spectrum than does the flat prior case. Note that the restriction on the upper limit for the 
mass parameters ensures relatively large production cross sections at the LHC for the case 
of the flat prior model sample. 

Once these points were generated, we demanded that they be consistent with a large 
number of both theoretical and experimental constraints in order to ensure that the model 
sets are valid to study. We mention the most important of these restrictions her^J (i) 
The spectrum is required to be tachyon free, color and charge breaking minima must be 
avoided, a bounded Higgs potential must be obtained and elect roweak symmetry breaking 
must be consistent. (ii) We impose a number of flavor and electroweak constraints arising 
from g — 2, b 57, B ru : B$ —> \i~ , meson-anti-meson mixing, the invisible width 
of the Z and Ap. (Hi) We demand that the LSP contribution to the dark matter density 
not exceed the upper bound determined by WMAP; note that the LSP is not required to 
saturate the measured relic density, leaving room for the existence of other dark matter 
candidates. Constraints from dark matter direct detection searches are also applied, (iv) 
We then include the restrictions imposed from the numerous direct searches for both the 
SUSY particles themselves as well as the extended SUSY Higgs sector at LEP. Here, some 
care was required as some of these searches needed to be re-evaluated in some detail due to 
particular SUSY model-dependent assumptions present in the analysis which we needed to 
remove, (v) Finally, the null results from a number of Tevatron searches are imposed. In 
addition to the Higgs searches, the most restrictive searches were found to be those hunting 
for stable charged particles [16] and those looking for an excess of multijet plus MET events 
fTT] . We note that in the latter case, the search strategies were designed for kinematics 
expected in mSUGRA-like models. We thus were forced to simulate them in some detail, at 
the level of fast Monte Carlo, for our full model set. At the end of this analysis chain, <~ 68.4/c 
models from the flat prior set survived this set of constraints, as well as a corresponding set 
of ~ 2.9k log prior models. These are the models that we will consider in our following 
analysis. 

We now turn our attention to the analysis procedure that we followed in generating 
and analyzing the signal events from sparticle production at the 7 TeV LHC. Throughout 
our analysis, we adhere to the search strategies developed by ATLAS p2] as closely as 
possible. In this reference, ATLAS considers 10 MET search channels, including selections 
where the minimum number of jets is 2, 3, or 4 and the number of leptons is 0, 1, or 2. In the 
dilepton case, opposite-sign (OSDL) and same-sign (SSDL) pairs are considered separately, 
with SSDL only being considered in association with 2 jets. Note that flavor tagged final 
states are not considered here. We consider 85 SUSY production processes that contribute 
to these 10 signatures. 

Accurate estimates of the SM backgrounds for the various channels are crucial to 
the validity of this study. We obtained details of the background distributions presented 
in Ref. [12] directly from the ATLAS SUSY Group [18]. These backgrounds were produced 

2 For full details, see Ref. [7] 
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with state-of-the-art Monte Carlo event generators and the full ATLAS detector simulation. 
Employing these ATLAS computed backgrounds in our analysis allows us to concentrate on 
generating and analyzing signal events for each of the ^ 71k parameter space points in our 
pMSSM model sample. 

2.1 Generation of the Signal Events 

The steps involved in the generation of the signal events are very similar to those detailed in 
Ref. [9j. Here we will briefly summarize the procedure and point out any differences in the 
present analysis. For the generation and analysis of events for a single model, the workflow 
is: 

1. The spectrum and decay table was generated with a modified [9] version of SUSY- 
HIT [T9]. 

2. The NLO cross sections for the 85 distinct SUSY production processes considered were 
computed using Prospino2.1 [20] and the CTEQ6.6M parton distribution functions [2T] 
(which were also used in the event generation). 

3. Using PYTHIA 6.418, events were generated, fragmented, showered, and hadronized 
for each of 85 SUSY production processes, with each process being weighted by its 
K- factor. 

4. Detector effects were simulated using an ATLAS-tuned version of the fast detector 
simulation PGS4 [22] with the default isolation cuts removed. 

5. The simulated events were then analyzed using the analysis cuts for the 10 ATLAS 
analyses listed above, as well as the isolation cuts described in Ref. [12]. 

We note that as in our previous work, a subset (about 1%) of the models suffered seri- 
ous enough errors that the Pythia event generation halted. These "PYSTOP" models are 
excluded from our results. 

2.2 Analysis Cuts 

For the reader's convenience, we provide here the full set of kinematic cuts for each analysis 
channel, summarizing the information given in Section 4 of Ref. [12j. All channels have a 
missing energy cut of E™ 1SS > 80 GeV, and all analyses except the SSDL channel have a 
transverse sphericity cut of St > 0.2. Table [T] summarizes the cuts for the searches with 
n = 2,3,4 jets which are independent of the choice of lepton channel except for SSDL, which 
will be described below. 
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Number of jets 


> 2 


> 3 


> 4 


Leading jet pt (GeV) 


> 180 


> 100 


> 100 


Other jets pt (GeV) 


> 50 


> 40 


> 40 


min. A0(jet i5 £;^ iss ) 


0.2,0.2 


0.2,0.2,0.2 


0.2,0.2,0 


E™ ss > fxM eS 


/ = 0.3 


/ = 0.25 


/ = 0.2 



Table 1: The kinematic cuts employed in the event selection: the cut on the pt of the leading 
jet, the pt of the other selected jets, the azimuthal angle between the selected jets and the 
missing transverse energy, and the missing energy as a fraction of the effective mass. 

We complete the description of the kinematic cuts by specifying the additional cuts 
that are specific to the various lepton channels. For the case with zero leptons, events are 
rejected that have at least one lepton with p T > 20 GeV. For the one-lepton channels, one 
lepton with pt > 20 GeV is required, no additional leptons with pt > 10 GeV are allowed, 
and the transverse mass Mt of the selected lepton and the missing energy vector must satisfy 
M T > 100 GeV. (The definition of M T can be found in Ref. [4J.) For the OSDL channel, 
exactly two leptons with p T > 10 GeV are required, and they must have opposite charge. 

As mentioned above, the SSDL channel has distinct cuts; unlike all the other analyses, 
the three different jet selection options specified in Table [j] are not employed. Instead, two 
jets with p T > 80 GeV are required. In addition, two leptons with p T > 20 GeV, same 
charge, and invariant mass mu > 5 GeV must be present, and there is a veto on additional 
leptons with p T > 10 GeV. The transverse mass of the leading lepton and the missing energy 
vector must satisfy M T > 80 GeV. 

Lastly, when performing our statistical analysis, the M e ^ cut is optimized for each 
channel, and for each pMSSM model, in steps of 400 GeV [12]. 

2.3 Statistical Procedure 

To compute the significance of the signal for each search channel, we follow the statistical 
procedure described in detail in Refs. [9j[4], which is that employed by the ATLAS collabo- 
ration. The probability that the expected background fluctuates to the number of observed 
events is computed assuming that the systematic error on the background is Gaussian and 
the statistical error is Poissonian. A significance of S > 5 is required for the observation 
of a signal. As mentioned above, the M e R cut is optimized for each channel, and for each 
pMSSM model, in steps of 400 GeV [12J. As was discussed in Ref. [9j, and will be further 
demonstrated here, the accuracy of the background estimation has a profound impact on 
the signal significance and the resulting search reach. In order to quantify this, we will 
present results assuming a 20, 50, and 100% systematic error on the background. We will 
also consider integrated luminosities of 0.1, 1, and 10 fb _1 . 
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Figure 1: The M e ff distribution for the 4 jet, 0(1) lepton analysis on the left(right) for the 
SU4 benchmark model. The green data points represent our analysis, while the blue line is 
the result from the ATLAS study [12]. The red shaded area represents the SM background. 



2.4 Comparison with ATLAS Benchmark Models 

The ATLAS SUSY group has published signal rates for a single mSUGRA benchmark point 
(SU4) in their study of Supersymmetry at the 7 TeV LHC [12]. It is imperative for us to 
check the results of our analysis against these published results for this benchmark model 
before proceeding to apply our analysis to the pMSSM model set. 

Here it is important to remind the reader that our SUSY signal generation, as de- 
scribed in detail in Ref. [9j, differs slightly from the procedure employed by ATLAS. In 
particular, the numerical programs used to compute the SUSY spectrum and decay tables, 
as well as event generation, are different. Furthermore, we use, by necessity, a fast detec- 
tor simulation as opposed to the ATLAS full GEANT-based simulation. Therefore, a small 
degree of discrepancy can be expected. The comparisons shown in Figures [T]{5j however, 
indicate that we are indeed able to faithfully reproduce the results obtained by ATLAS for 
this benchmark model for all of the various inclusive analyses. 
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Figure 2: The same as Figure [TJ except for the 3 jet, 0(1) lepton analysis on the left(right). 
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Figure 3: The same as Figure [TJ except for the 2 jet, 0(1) lepton analysis on the left(right). 
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Figure 4: The same as Figure [TJ except for the 4(3) jet, OSDL lepton analysis on the 
left (right). 
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?ure 5: The same as Figure[IJ except for the 2 jet, OSDL(SSDL) analysis on the left (rig 
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3 Results of the 7 TeV Analysis 



In this Section, we relate the results of our study on the effectiveness of the ATLAS 7 TeV 
E™ lss analyses in detecting our pMSSM model sample. We first discuss the impact of the 
size of the background systematic errors on SUSY searches, and then turn to the discovery 
coverage of the pMSSM. We examine the characteristics which cause some models to be 
undetectable as well as study the effects of modifying the ATLAS SUSY analysis cuts. We 
remind the reader that our sample of ~ 70k models is not intended to be a full description of 
the 19-dimensional pMSSM parameter space. However, the sample does contain numerous 
models which exhibit properties that are quite different than those expected in mSUGRA 
and thus provides insight into general features of the full MSSM. 

3.1 Influence of Background Systematic Errors 

As mentioned above, the size of the SM background systematic errors plays an important role 
in the ability of the ATLAS E™ 1SS searches to discover Supersymmetry, including the pMSSM. 
This is not surprising as the number of signal events necessary to reach S = 5 critically 
depends upon both the size of the estimated background itself as well as the background 
uncertainty. For a fixed systematic uncertainty, search channels with large backgrounds 
clearly require a large number of signal events in order to claim a discovery. To get a feel 
for this in the case of the ATLAS E™ 1SS analyses studied here, we determine the necessary 
number of signal events to reach the S = 5 level in each analysis as function of the fractional 
background uncertainty. We remind the reader that the SM backgrounds for each channel 
were supplied to us by the ATLAS SUSY working group (THJ • In performing these calculations 
we exactly follow the discussion as given by ATLAS in Ref . [4J . Our results are displayed in 
Figs. [6] and [7] for the ten ATLAS E™ 1SS channels assuming 1 fb _1 of integrated luminosity. 
Here, we see the number of signal events that are required to obtain the discovery criterion 
of S = 5 for various values of the final M e ff cut. In the case of the njOl channel, which has 
the largest SM background, we note that the required number of signal events is particularly 
large and is quite sensitive to the value of the M e ^ cut. Note that as the systematic error 
increases, the number of required signal events can rise drastically, in some cases by an order 
of magnitude or more. In particular, the difference between a reasonably low 20% systematic 
error and taking a 0% error (i.e., ignoring this effect) is substantial and theoretical analyses 
that do not include this error are thus wildly optimistic. 

We will use these numerical results in our subsequent analyses of the pMSSM model 
coverage in these £7™ lss -based searches in the next subsection. They indicate the importance 
of reducing background systematic errors in order to increase the coverage of new physics 
parameter spaces. 
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N s required to get 5o discovery with various M eff cuts for njOl for 1 fb 




1 1 1 1 1 1 1 

20 40 60 80 100 

Systematic Error 

Figure 6: Number of events required to reach the S = 5 level of discovery as a function of the 
fractional systematic error in the SM background for the ATLAS E™ 1SS searches for various 
values of the M e ff cut. The results for the njOl and njll searches are shown in the top and 
bottom panels, respectively. The curves are color coded from top to bottom as indicated 
in the legend. For higher values of the M e ff cut, we see that the curves are essentially 
indistinguishable, lying on top of one another. 
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N s required to get 5o discovery with various M eff cuts for njosdl for 1 fb 
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N s required to get 5o discovery with various M e ff cuts for 2jssdl for 1 fb" 
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Figure 7: Same as the previous figure but now for the njSSDL(top) and 2jSSDL (bottom) 
search channels. 
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3.2 pMSSM Model Coverage 



We now run each of our pMSSM models through the analysis chain described above. The 
first question we address is how well do the various search analyses cover the pMSSM model 
sample, or, more precisely, what fraction of these models can be discovered (or not) by 
these searches. Further, we also determine which of the analyses provide the best model 
discovery capabilities. Clearly the answers to these questions will be highly sensitive to 
the assumed values of both the integrated luminosity and the estimated SM background 
uncertainty. Figure [8] shows the fraction of the pMSSM models that can be discovered with 
S > 5 in each of the ATLAS E™ 1SS channels for the flat prior model set as a function of the 
integrated luminosity assuming three different choices for the background systematic error. 
The corresponding results obtained in the case of the log prior model set can be found in 
Fig. [9| Here, we again emphasize that our pMSSM sample is not meant to provide full 
coverage of the 19-dimensional parameter space (such coverage would be computationally 
prohibited). However both the very large number of models in our pMSSM sample, and 
the distinct characteristics they possess, make this sample an ideal test bed for this set of 
mSUGRA designed search strategies. Our fractional results based on our pMSSM model set 
are thus indicative of the behavior of the MSSM under these search routines. 

These figures reveal a number of interesting results: (i) The size of the background 
systematic errors makes a significant impact on model coverage for all search channels and 
integrated luminosities. For the searches with significant SM backgrounds, i.e., the njOl and 
njll channels, variation in the background uncertainty leads to substantial changes in the 
capability to observe the model sample. The search least affected by systematics is 2jSSDL 
since the backgrounds in this case are quite small. This behavior confirms the results of the 
previous subsection, (ii) The model coverage in almost all cases is significantly better for 
the flat prior model set than for the log prior sample. This, too, is not very surprising since 
the masses of the sparticles in the log prior case extend out to much larger values and the 
sparticle spectrum is generally more compressed in this set [7J. The latter leads to softer jets 
and leptons in the corresponding cascade decays which have a more difficult time passing 
the analysis cuts. (Hi) For all values of the background systematic error, the njOl channels 
yield the best model space coverage with 4j01 affording the best discovery opportunity. In 
fact, we see that the channels which require more jets to be present have a better chance 
of being observed for the njOl and njll searches. For the flat prior sample, the 4j01 analysis 
with low background systematics is observed to cover a very large fraction of the model 
set by itself once significant integrated luminosities are obtained, (iv) As the number of 
leptons required to be present in the final state increases, the model coverage is found to 
decrease significantly, especially for smaller values of the integrated luminosity. This is due 
to the fact that the branching fractions for leptons to appear in squark and gluino induced 
cascade decays are generally not very large in our model sample, as we have seen in our 
earlier work [9]. (v) Independently of the specific E^ 1SS search, as the background systematic 
errors become large, the pMSSM model coverage is seen to increase more slowly with the 
integrated luminosity. Some of the search channels nearly saturate at high luminosity due 
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Figure 8: Fraction of flat prior pMSSM model set that can be observed with S > 5 in the 
njOl(top-left), njll(top-right), njOSDL(bottom-left), and 2jSSDL(bottom) search channels 
as a function of the integrated luminosity. The solid(dashed, dotted) curves in each case 
correspond to n=4(3,2), respectively for the njOl, njll, njOSDL channels. The red(green, 
blue) curves correspond to background systematic uncertainties of 20(50, 100)%, respectively. 
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Figure 9: Fraction of log prior pMSSM model set that can be observed with S > 5 in 
the njOl(top-left), njll(top-right), njOSDL (bottom-left), and 2jSSDL(bottom-right) search 
channels as a function of the integrated luminosity. The solid(dashed, dotted) curves in 
each case correspond to n=4(3,2), respectively for the njOl, njll, njOSDL channels. The 
red(green, blue) curves correspond to background systematic uncertainties of 20(50, 100)%, 
respectively. 
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to the large background uncertainties and thus become systematics dominated. 

As discussed above, the final step in the ATLAS E™ 1SS analyses is to apply a cut 
on M e ff, where the particular value of M e ^ that is chosen (in units of 400 GeV) is the 
one that maximizes the signal significance given the SM background and its corresponding 
uncertainty. This choice not only depends upon the particular channel but also on the 
amount of integrated luminosity. It is important to note that if this cut is taken to be too 
large when maximizing the signal, then the analysis will be very sensitive to the detailed 
shape in the tails of both the signal and expected background distributions, especially with 
higher luminosities. This happens when there are very few events with large values of M e ff, 
e.g., 2 TeV. In this situation, small fluctuations in the SM background and/or SUSY signal 
expectations due to limited Monte Carlo statistics can lead to inconsistencies in whether a 
given model is observable in a specific analysis or even whether or not it would be detected 
overall. This effect only occurs in the case of the search analyses designed for the 7 TeV run, 
as the M e ff cut applied in the planned 14 TeV ATLAS analyses was fixed at relatively low 
values. 



Figure [10] shows the optimized value for the M e ^ cut for the njOl analyses, as an 
example, for the flat prior model sample. The three things we see here are: (i) for large 
background systematic errors, a harder M e ff cut is required to optimize the search, (ii) as the 
number of required jets in the final state decreases, the strength of the cut can be reduced. 
Both of these results are also found to hold for the njll and njOSDL searches although the 
M e ff cut itself turns out to be less important as the number of required leptons in the final 
state increases. (Hi) Given the warning about distribution tails in the discussion above, it 
is a welcome result to see that in the majority of cases only a moderately strong M e ff cut 
is required to optimize the signal significance. Note that there are some cases where the 
M e ff cut does not contribute very much to increase the significance of the SUSY signal; this 
happens in particular for scenarios where the background is low and the effect of systematic 
uncertainties is not very significant. 

It is interesting to evaluate the fraction of models that can be discovered in multiple 
analyses. This is important to consider as, first, it is valuable to verify the discovery of new 
physics in more than one channel, and second, the availability of multiple discovery channels 
admits for the possibility of further studies that will allow for, e.g., the extraction of sparticle 
masses. To be specific, Tables [2j [3] and [4] show the fraction of pMSSM models which are 
observed in (exactly) n channels assuming a SM background systematic error of 20, 50, and 
100%, respectively. Note that the distribution shifts towards more models being observed 
in multiple channels as the integrated luminosity increases and the background systematic 
error decreases, as expected.. 

By combining our results for these E^ 1SS searches, we can also determine the fraction 
of pMSSM models that are undetected in all of the 7 TeV search analyses designed by 



ATLAS; this corresponds to the case of n = in these Tables. Figure 11 presents the 



fraction of pMSSM models which are undetected in all of the search channels as a function 
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Histogram of the best M eff cuts for 4jOI for 1 fb 1 
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Histogram of the best M eff cuts for 3jOI for 1 fb 1 



Histogram of the best M eff cuts for 2jOI for 1 fb 1 





800 1200 1600 2000 
Meff cuts (GeV) 



2400 2800 



800 1200 1600 2000 
Meff cuts (GeV) 



2400 2800 



Figure 10: Optimized M e ff cut for the flat prior model set for the 4j01(top), 3j01(bottom left) 
and 2j01 analyses (bottom right). The red(green, blue) histograms correspond to background 
systematic uncertainties of 20(50, 100)%, respectively. An integrated luminosity of 1 fb -1 
has been assumed in these figures for purposes of demonstration. 
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Table 2: The fraction of models that are observed in (exactly) n E™ 1SS search channels 
assuming a SM background systematic error of 20%. 
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Table 3: Same as the previous Table but now assuming a SM background systematic error 
of 50%. 
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Table 4: Same as the previous Table but now assuming a SM background systematic error 
of 100%. 



of integrated luminosity for both the flat and log prior sets. In the flat prior case we see 
that as the integrated luminosity increases, the model coverage substantially improves, and 
approaches (or exceeds) ~ 95% for 10 fb _1 with SM background uncertainties of 50% or 
less. In the log prior case, the improvement in pMSSM model coverage as the luminosity 
increases is much more gradual as we expected. However, even in this case, at high integrated 
luminosities substantial model coverage is seen to be obtainable at 7 TeV. 

This figure also shows the important playoff between increasing the integrated lu- 
minosity and decreasing the SM background systematic error in terms of pMSSM model 
coverage. (Of course, increased luminosity often results in decreased systematic errors, up 
to a point.) For example, it is interesting to compare the effectiveness of the analyses for 



the flat prior set assuming SB = 100% and C = 5(10) fb 1 with other values. Figure 11 
shows that taking SB = 50% and C = 0.65(1.4) fb" 1 or SB = 20% with C = 0.20(0.39) fb" 1 
produces essentially identical model coverage. This demonstrates that small reductions in 
the SM background uncertainty can be worth a significant amount of increased integrated 
luminosity in terms of pMSSM model coverage. 

Another very important message to take home from this figure is that for the inte- 
grated luminosity collected during the 2010 LHC run at 7 TeV (~ 35 pb _1 for ATLAS), 
some significant fraction of these pMSSM models should already have been observed at the 
S = 5 level. Explicitly, for a background systematic error of 20(50,100)% we find that 
46(28,15)% of the flat prior model set should have been discovered; in the log prior case 
the corresponding results are found to be 30(24,14)%. Since these are discovery results, an 
even greater portion of the pMSSM model sample would be expected to be excluded by these 
analyses. This shows the incredible power of going from Tevatron energies up to the 7 TeV 
LHC in performing searches for massive objects, such as SUSY sparticles, even when only 
small amount of integrated luminosity is available. 



19 



to 



CD 



n 

cd 
o 

CD 



100 




0.05 0.10 0.20 0.50 1.00 2.00 

-1> 



5.00 10.00 



Luminosity (fb ) 





3 

i 

•j: 
■ — i 



CD 



loo n — i i i i 



-i — i — i — i i i 1 1 



-i — i — i — i i i 




20 — 



_i i i i i i i i 



_i i i ' ' ' ■ 



0.05 0.10 0.20 0.50 1.00 2.00 5.00 10.00 

Luminosity (fb -1 ) 



Figure 11: Fraction of flat prior (top) or log prior (bottom) pMSSM model sets which are 
undetected after combining all of the ATLAS E™ 1SS search analyses. The red(green, blue) 
curves correspond to background systematic uncertainties of 20(50, 100)%, respectively. 
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3.3 Why are Models not Detected by the ATLAS Ef ss Searches? 

Here, we investigate the main reasons why some pMSSM models are not discoverable in the 
ATLAS E™ ss searches at 7 TeV. We addressed this question in some detail in our earlier 
work for the case of the analyses designed for the 14 TeV LHC [9j, so our discussion here 
will not be as extensive. Clearly, many of our previous results will carry over qualitatively 
into the present 7 TeV analysis. 

There are multiple explanations as to why some pMSSM models are undetected by 
the £j? lss searches, the most obvious one being small production cross sections for the colored 
sparticles that initiate the familiar SUSY cascades. As an example, we note that for our 
pMSSM models in the flat prior set, the cross sections for the production of gluino and 
squark pairs are found to cover an enormous range of several orders of magnitude as can 
be seen in Fig. [12] (recall that the upper limit on sparticle masses in our flat model set is 
~ 1 TeV). Here we see that the large (or small) values of the gluino pair cross section is 
completely uncorrelated with the corresponding values for first generation squarks within 
a particular model. Furthermore, by summing over all of the QCD production channels 
involving gluinos and/or first generation squarks (i.e., gg, gq ) qq and qq*) we obtain an 
approximate handle on the total overall rate for SUSY production which we see ranges over 
four orders of magnitude. Note that for any particular value of the squark or gluino mass, 
the corresponding production cross section itself can vary by up to an order of magnitude 
or more depending upon the remainder of the pMSSM model spectrum. 

While it is certainly clear from this figure that some models have too small a cross 
section to be discovered, the bottom-right panel indicates that this cannot be the entire 
explanation. Here we show the search significance, S, of the 4j01 channel (as it is the most 
powerful channel in terms of discovery capability) as a function of the total NLO QCD 
production cross section assuming C = 1 fb _1 and 5B = 50%. Here we observe that (i) 
there are models with cross sections ~ 20 pb which are missed by this analysis, while (ii) 
there are models with cross sections ^ 100 fb which are discovered. (Hi) For any given 
value of the cross section, the range of the significance is large and can be up to two orders 
of magnitude or more. This validates the claim that there are reasons other than small 
production cross sections that render models unobservable by these E^ 1SS analyses, (iv) 
For any given gluino mass there is a strong correlation of the signal significance with the 
mass of the LSP. Clearly when these two masses are close the average pt of the jets will 
be softer and this will make it more difficult to pass analysis cuts. Also if the LSP mass is 
large then that implies even larger squark and gluino masses that will result on average in 
smaller production cross sections. In the case of the njOl analyses, a larger number of signal 
events is required for discovery due to the sizeable SM backgrounds and hence such models 
will be missed by these analyses. Visibility will then require the production of leptons with 
significant pt in cascade decays in order to pass the lower background njll, njOSDL and 
2jSSDL searches. Unfortunately, lepton branching fractions are low in these cascades in our 
model sample (since, e.g., sleptons are heavy) and thus some models will be missed entirely. 
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Figure 12: (Top left ) NLO first generation squark and gluino pair production cross sections 
at y/s = 7 TeV as a function of their masses for the flat prior model set. The green(red) 
points represent the gluino (squark) cross sections. (Top right) Correlation of the squark pair 
and gluino pair NLO cross sections in the flat prior set. Each point represents one model. 
(Bottom left) Total NLO QCD production cross section distribution for the flat model set. 
(Bottom right) Search significance of the 4j01 analysis as a function of the total NLO QCD 
production cross section assuming C — 1 fb _1 and 8B = 50%. The solid line highlights the 
S = 5 discovery level and the color code reflects the mass of the LSP. 
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Figure 13: Significance of the 4j01 search for the flat prior model set as a function of the 
gluino(left) and lightest squark(right) masses. C = 1 fb _1 and 5B = 50% have been assumed. 

Of course, models that are lepton-rich will automatically fail all the njOl analyses since they 
veto events with high p T leptons, but will be picked up by the searches containing leptons. 

Though the production cross section for SUSY particles is reasonably well correlated 
with their masses, we can ask more directly if larger sparticle masses lead to their non- 
observation in these searches. Figure p~3| shows that, indeed, models with lighter squarks or 
gluinos tend to lead to signals with greater significance in the 4j01 channel. However, as we 
can also see from this Figure, this is not true universally, e.g., there are many models with 
gluino (lightest squark) masses below 300 (200) GeV that have S < 5, while conversely there 
are models with 1 TeV gluinos that have S > 5. We see that for any given squark or gluino 
mass the value of S can vary significantly. The top panel of Fig. [14] displays this property 
even more strongly where we see that these results hold even when all of the E^ 1SS searches 
are combined. This Figure shows the set of flat prior models that are unobservable in all 
the search channels in average light squark mass-gluino mass plane. Note that there are 
a significant number of these models which contain light squarks and gluinos. Thus while 
the masses of the colored sparticles do play an important role in model observability clearly 
there are additional important factors. 



In the lower panels of Fig. 14 we see that the mass splitting between squarks and/or 
gluinos and the LSP can play an important role in determining model observability as was 
first noted in Ref. [6j and was seen explicitly in our earlier work on the generation of the 
pMSSM models [7J and the 14 TeV ATLAS SUSY analyses [9j. The obvious reasoning here 
is that as the degeneracy in the spectrum increases and mass splittings become smaller, the 
values of, e.g., the pt of the jets, will be reduced so that it will be more difficult to satisfy any 
of the analysis cuts. These figures show this result explicitly. Note, however, that, e.g., in the 
case of light gluinos with small gluino-LSP splittings there are many models which are still 
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discoverable in the 4j01 channel. The reason for this is that while the efficiency^ for passing 
the 4j01 analysis cuts may be quite low for small mass splittings, as seen in Fig. [15J the cross 
section to produce the lighter gluinos/squarks is very large and more than compensates for 
these low efficiencies, especially if there are any additional hard jets in the event from ISR. 
However, there are numerous unobservable models that have larger raw sparticle production 
cross sections than observable models with somewhat similar spectra; the difference then 
being in their respective abilities to pass the necessary analysis cuts. 

Another cause for models being undetected is the occurrence of detector-stable spar- 
ticles at the end of gluino or squark induced decay chains instead of the LSP [2J. This 
happens with reasonable frequency in both the 7 and 14 TeV analyses. In such cases, the 
amount of E™ 1SS that is produced is substantially decreased which reduces the capability of 
the relevant models to pass any E™ 1SS analysis requirements. Most commonly, these sparti- 
cles are actually long-lived charginos that are reasonably degenerate with the LSP in wino- 
or Higgsino-like LSP scenarios. In these cases, searches for long-lived sparticles, as discussed 
below, will be a important supplement to the conventional E™ 1SS searches. Of course, a loss 
of the E™ ms signature can happen in other ways. For example, if the initial squark or gluino 
produces a very long decay chain then the particles produced at the end of such a chain 
will be somewhat soft. In some cases this may lead to the inability to pass the necessary p T 
and/or E^ 1SS requirements for the various £^? lss -based searches and the model will not be 
observed. Such long decay chains were shown to occur with a reasonable frequency in our 
earlier work [9j and can be a contributor to models failing to pass the analyses requirements. 

There are, of course, other reasons that prevent models from being discovered. As 
noted above, subtleties in any sector of the sparticle spectrum can make a significant differ- 
ence as to whether a given model is observed by various analyses. Here, we will discuss a 
couple of examples where this occurs. As in Ref. [9], the approach we follow is to compare 
a model which fails to be observed in all search channels to one with a similar spectrum 
(dubbed a 'sister' model) that is detected in at least one channel and then examine the dif- 
ference between them. For this study, we concentrate on the more difficult cases by taking 
the flat prior model sample and assume C = 10 fb _1 and 8B = 20%; this leaves only ~ 670 
models that are not observed. In order to avoid the statistical issues associated with the 
tails of the M e ^ distribution discussed above, and to further reduce this model set to a more 
manageable size, we will only consider models whose optimized M e ^ value is < 800 GeV. 
We note that a large fraction of models in this set have relatively heavy LSPs with masses 
in excess of 400 GeV. We now discuss two brief examples of these comparisons. 



Figure 16 shows a comparison of two similar models, 8944 (observed in the 3,4jOSDL 
channels) and 21089 (missed by all analyses), that have comparable total colored sparticle 
production rates (3.4 and 4.6 pb, respectively). Both of these models are not observed in 
the njOl searches since the lighter squarks are too close in mass to the LSP to produce hard 
jets. The gaugino sectors of these two models are quite similar (with the LSP and \2 being 



3 Here, efficiency is defined as the fraction of generated signal events that pass the analysis cuts. 
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Figure 14: (Top) Flat prior models that are unobservable in all of the £ , ™ lss -based search 
analyses in the average light squark -gluino mass plane. (Bottom) Flat prior models that 
pass(green) or fail(red) the 4j01 analysis in the gluino (lightest squark) mass vs gluino- 
LSP (lightest squark-LSP) mass splitting plane in the left (right) panel. C — 1 fb -1 and 
5B = 50% have again been assumed. 



25 



efficiency us mg-L£P for pass models 




400 500 
mg-LSP or sq-L£P 



gure 15: The efficiency for passing the 4j01 analysis cuts as a function of the (green) gluino- 
(red)squark-LSP mass splitting. Here, efficiency is defined as the fraction of generated 
pial events that pass the analysis cuts. 
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Figure 16: The top panels compare the leading jet (left) and leading lepton (right) pt spectra 
for models 8944 and 21089. The bottom panels show the sparticle spectra for these two 
models. 
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Figure 17: The top panels compare the leading(left) and secondary(right) lepton p T spectra 
for models 9781 and 20875. The bottom panels show the sparticle spectra for these two 
models. 



Higgsino-like and X3 being bino-like), while their colored sparticle spectra are somewhat 
different. In either model the decay of X3 allows for the OSDL production through an 
intermediate slepton which has sufficient pt to pass the analysis requirements. However, 
while 8944 has a ur with mass below that of the gluino (which is light enough to give a 
reasonable cross section), allowing for the decay into X35 on ly dR is (sufficiently) heavier 
than the X3 i n model 21089. In this case, is also much more massive than the gluino, 
through which it will dominantly decay, and so it will not have a large enough branching 
fraction into Xs to produce the OSDL signature. 

Figure [Tt] compares models 9781 (discovered in the 2jSSDL channel) and 20875 (com- 
pletely missed) that have total colored sparticle production cross sections of 1.3 and 1.1 pb, 
respectively, but yet produce too few hard jet plus E^ 1SS events to be found in the njOl 
channels due to spectrum compression. Model 9781 is quite interesting as the charginos and 
lightest three neutralinos are highly mixed combinations of winos and Higgsinos. In this 
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model, ur (which is relatively light) decays to j + X2 with a ~ 98% branching fraction. Since 
X2 has a large bino content, it decays ~ 95% of the time through sleptons which subse- 
quently decay directly to the LSP with a branching fraction of ~ 45%. Thus model 9781 can 
easily populate the leptonic final state and since the neutralinos are Majorana fermions, the 
2jSSDL final state becomes accessible, a channel with an extremely small background. On 
the other hand, model 20875 does not allow for the generation of a typical leptonic signal. 
This is because the e, jl are quite heavy so that the neutralinos only allow for decay to r 
leptons via an intermediate on-shell f . This is correlated with the relative lightness of the 
f as well as the Higgsino-like LSP. In addition, the lightest squarks decay primarily directly 
to the LSP which will not produce any high-p T leptons. 

These two examples demonstrate that the full sparticle spectrum may conspire to 
render a model undiscoverable in the £j? lss -based analysis suite, even if the colored sparticle 
cross section is large. Discovery is not based on the value of the squark and/or gluino masses 
alone and blanket limits that claim Tfiq g are ruled out below some value cannot be set. 



3.4 Detector Stable Sparticles in Cascades 

As we mentioned above, one of the reasons that some pMSSM models may not be observed 
in the ATLAS E^ lss analyses is that squark and gluino cascade decays can sometimes lead 
to a final state with low E™ 1SS - In many cases this is due to the existence of long-lived 
sparticles, usually charginos, which appear with sizable branching fractions in such cascade 
chains and are essentially detector stable. Particles that decay outside the detector when 
produced at the LHC would provide a dramatic signal of new physics (see [23] and references 
therein). In fact, data from the LHC are already extending the mass limits on such detector- 
stable sparticles [24] . Since the inclusive E^ 1SS analyses discussed above do not consider such 
sparticles, we will briefly sketch the 7 TeV discovery prospects for detector-stable sparticles 
in our model sets. As discussed in our earlier work [9], the existence of such long-lived states 
is relatively common in our pMSSM model sample as can be seen in Fig l8|^ 



Our estimation of the 7 TeV LHC mass reach for each long-lived sparticle is shown 
in Table [5] for the specified integrated luminosities, assuming direct pair production. These 
results are deduced from Figure 1 of [25] by taking the geometric mean of the search reach 
for a 5 and 10 TeV LHC at the specified luminosities and are interpolated to the lumi- 
nosities considered here where necessary. These results are somewhat conservative as only 
detector-stable sparticle production in the hard process is considered; additional detector- 
stable sparticles could be produced through cascade decays as discussed below. 

The number of detector-stable sparticles of various species in our model sample is 



4 Recall that we will define a particle to be detector-stable if its unboosted decay length in at least 20m. 
Note that typical values of 7/3 for long-lived particles resulting from cascade decays are in the 2-3 range. 
Further note that the dependence of the number of detector-stable sparticles of various species in this model 
set on the value of the decay width r sta bie that is assumed is discussed in detail in [9]. 



29 



Sparticle 


LHC Reach 100 pb" 1 GeV 


LHC Reach 1 fb" 1 


LHC Reach 10 fb" 1 


X + (Wino-like) 


206 GeV 


264 GeV 


334 GeV 


X + (Higgsino-like) 


153 GeV 


204 GeV 


267 GeV 


f 


79 GeV 


109 GeV 


146 GeV 


t 


294 GeV 


363 GeV 


441 GeV 


9 


563 GeV 


654 GeV 


751 GeV 



Table 5: Approximate 7 TeV LHC search reaches for detector-stable sparticles of the given 
species with lOOpb" 1 , 1 fb" 1 , and 10 fb" 1 [25!. 
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Table 6: The second column from the left gives the number of detector-stable sparticles of 
various types in our model set. The next two columns show the number of such sparticles that 
will not be discovered after 100 pb" 1 , 1 fb -1 , and 10 fb -1 at the 7 TeV LHC, following [25]. 



shown in Table [6j This Table also shows the number of detector-stable sparticles which 
will not be discovered at LHC with 100 pb -1 and 1 fb -1 of integrated luminosity, using the 
approximate mass reaches presented in Table [5j We assume here that the mass reach is 
roughly generation-independent and that it is the same for, e.g., stops and sbottoms. This 
assumption is reasonable except where there could be significant i-channel production for 
the first or (to a lesser extent) second generation, for instance in the case of up or down 
squarks. 

We now focus on the specific case of long-lived charginos, which are by far the most 
common long-lived sparticles in our model sets. If the production cross section for colored 
sparticles (times relevant branching fractions into charginos) are sufficiently large, these sta- 
ble charginos should be found in searches for (effectively) stable charged particles occurring 



at the end of a cascade decay chain. In Fig. 18 we display the estimated value for aB for the 



production of detector-stable charginos in cascade decays in our flat prior set. (Note that 
this does not include the direct contribution arising from direct chargino pair production.) 
Here, we see that roughly ~ 84% of models with detector-stable charginos lead to aB values 
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in excess of 10 fb at 7 TeV and so we expect them to be observable in the upcoming run 
of the LHC. In this estimation, we assumed the largest contribution to the production cross 
section arises from the production of gluino and light squark species. Using the decay tables 
generated for each model with detector-stable charginos, we calculated the branching frac- 
tion for the gluino and light squarks to produce a stable chargino at the end of each possible 
decay chain and then weighted them by their corresponding production cross-sections. Note 
that mass information for neither the mother particles or the daughter chargino is used to 
indicate how likely it is that the chargino will pass the trigger criteria for detection. As we 
used PYTHIA to compute the LO cross-section for the light squarks, we do not separately 
generate the production for the various light squark species. Thus in order to obtain our 
estimate, we make the assumption that the overall cross section is 100% dominated by that 
arising from the lightest squark. This assumption will break down when the squarks are 
nearly degenerate, especially when their branching fractions to charginos are vastly different 
due to the complexities in the gaugino sector. Many of our models have a large production 
cross-section for colored sparticles, but the relevant branching fractions to charginos can be 
simultaneously quite small. This can result in very small overall production rates for stable 
charginos well below ~ 10 fb. 

3.5 SUSY Mass Scale From M eff 

In our earlier work on the 14 TeV ATLAS E™ 1SS analyses, we demonstrated that the rela- 
tionship between M e ff and the mass of the lightest colored sparticle found in mSUGRA, i.e., 
M e ff ~ 1.5mLCP (where LCP stands for Lightest Colored Particle), proposed long ago [26] 
does not necessarily hold in the pMSSM. This possible relationship is important as it might 
be used to get the first handle on the overall mass scale of the sparticle spectrum. Here, we 
briefly note that this result remains valid for the 7 TeV ATLAS E™ 1SS analyses as can be 
seen in Fig. [l9j For both the 4j01 and 2j01 channels we see explicitly that the values of M e ff 
lie mostly above the expected value of 1.5mLCPj especially in the low sparticle mass region. 
However, for lightest colored sparticle masses in excess of ~ 550 — 600 GeV we see that, 
indeed, the relationship M e ^ ~ 1.5tolcp provides a fairly good estimate in both of these 
E™ 1SS searches. 

3.6 Modifying ATLAS SUSY Analysis Cuts 

Given the properties of the various sparticles in our model sets, we can try to determine 
whether the canonical cuts employed in the ATLAS i?^ lss -based search analyses can be 
strengthened to reduce SM backgrounds without any significant loss in the coverage of our 
pMSSM model space. This is certainly a non-trivial issue and the structure of our analysis, 
being based on the fixed ATLAS E™ 1SS analyses cuts, is not directly set up to obtain com- 
pletely definitive answers. However, it is possible to make some reasonable estimates based 
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Figure 18: (Top) Number of models from both the log and flat prior sets combined having a 
charged sparticle with an unboosted decay length above a given value. (Bottom) Distribution 
of the estimated cross section times branching fraction for the production of detector-stable 
charginos in cascade decays in the flat prior model set. 
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M eff vs. Minimum colored sparticle mass for 4jOI 



M eff vs. Minimum colored sparticle mass for 2jOI 



1600 
1400 
1200 
1000 
800 
600 
400 
200 





Fail 

Pass 



100 200 300 400 500 600 700 800 900 1000 1100 
Minimum colored sparticle mass (GeV) 



1600 r- 

1400 - 

1200 - 

1000 - 

800 - 

600 - 

400 - 

200 - 




Fail 



100 200 300 400 500 600 700 800 900 1000 1100 
Minimum colored sparticle mass (GeV) 



Figure 19: Correlation between the value of M e ff and the mass of the lightest colored spar- 
ticle for the 7 TeV 4j01(left) and 2j01(right) ATLAS Ef ss channels. The red(green) points 
correspond to flat prior models which are missed (found) in these two search analyses. 



on the information that we do have available. We will concentrate on the three njOl analyses 
as they generally provide the greatest pMSSM model coverage and have large statistics. The 
most important kinematic quantities for these searches are the requirements on the leading 
jet pt and the required amount of E^ 1SS . Here, we make use of the average values of the 
distributions in these quantities for our pMSSM model sample, as well as the correspond- 
ing fitted width of the part of the distribution below this average value for pre-selected 
events. This information then provides us with an estimate of where these two kinematic 
distributions 'turn on,' which we take to be the average value minus this width, on the low 
energy /momentum side below their peak average values. 

First consider the cut on the leading jet, , for the 2(3,4)j01 analyses; ATLAS chooses 



the value for this cut to be 180(100,100) GeV, respectively. Fig. 20 shows the distribution 
of the 'turn-on' pt values for these three ATLAS analyses obtained from analyzing the flat 
prior model set. For the 2j01 analysis, we see that the lower edge of the 'turn-on' values 
lies somewhat below the cut value of 180 GeV from which we can conclude that this cut is 
already reasonably hard and cannot be increased without a loss of model coverage. However, 
for both the 3j01 and 4j01 analyses we instead observe that the 'turn-on' values lie above those 
of the ATLAS cuts by ~ 20 GeV suggesting that that the p? x cut in these two channels may 
be raised without impacting coverage rates for the pMSSM. 

Similarly, the nominal E^ 1SS cut imposed by ATLAS for the njOl analyses is 80 GeV. 
However, there is an additional subsequent cut imposed by ATLAS based on the value of 



M eff , i.e., Ef ss > /M eff where / = 0.30(0.25,0.20) for the n = 2(3,4) analyses. Fig. [21 
shows the distributions for the 'turn-on' values of E™ 1SS for these three channels employing 
the flat prior model set. Here we see that the lower edge of this distribution occurs at 
~ 160(130, 120) GeV for n = 2(3,4), respectively. This is suggestive that the nominal E™ 1SS 
cut made by ATLAS may be increased for these three analyses without losing significant 
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Figure 20: Approximate location of the lower edge of the leading jet pt distribution as 
a function of the number of preselected events employing the flat prior model set for the 
2j01(top), 3j01(lower left) and 4j01(lower right) ATLAS search analyses as discussed in the 
text. The color code reflects the number of models. 
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Figure 21: Same as in the previous figure but now for the E^ 1SS distribution. 



pMSSM model coverage. 



4 Implications of the 7 TeV Run 



In this section, we explore some implications of a null search for Supersymmetry at the 7 
TeV LHC. We examine the degree of fine-tuning that would be placed on our pMSSM model 
sample and we discuss the resulting expectations for sparticle production at a 500 GeV and 
1 TeV Linear Collider. 



4.1 Fine-tuning in the Undiscovered pMSSM Models 

As has recently been discussed in the mSUGRA/CMSSM context [271 EH], it is apparent 
that if SUSY signatures are not discovered at the 7 TeV LHC as the integrated luminosity 
accumulates it is likely that the SUSY parameter space must become more fine-tuned, and 
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hence more problematic as a solution to the hierarchy problem. Since we know which models 
in our sample are discoverable (or not) by the ATLAS E^ 1SS search analyses, we can ask 
whether this same result also holds in the case of our pMSSM model setsj^] 

Figure [22] display the results of this analysis assuming a background systematic error 
of 50% for both the flat and log prior model samples. Clearly, in the flat prior case, one sees 
that as the integrated luminosity increases and more models can be discovered by ATLAS, 
those remaining yet undiscovered tend to be more fine-tuned as expected. In other words, 
the fractional loss of models from the full distribution occurs more rapidly with increasing 
luminosity for models with smaller amounts of fine-tuning. This is not too surprising as, 
overall, models with less tuning tend to have lighter SUSY sparticle spectra and are thus 
more easily discovered at the LHC. On the other hand, the results from the log prior model 
set appear to be affected somewhat differently in that the overall shape of the fine-tuning 
distribution does not appear to change very much by removing models that should have been 
already been discovered by ATLAS as the luminosity increases. In this case, we see that there 
is not much of an increase in the amount of fine-tuning as the set of undiscovered models 
shrinks. This represents one of the few apparent differences between these two different 
model sets. This can be explained by the fact that while the log prior model set tend to 
have light sparticle spectra (though they extend out to larger mass values than do those 
for the flat prior models) and are thus less fine-tuned to begin with, these same spectra 
are generally compressed making these models more difficult to discover at the LHC as was 
discussed above. This would imply that the models missed by the ATLAS E™ 1SS analyses 
in the log prior case are generally not much more fine-tuned than those appearing in the 
originally generated model set. Thus we find that the amount of fine-tuning that remains in 
the LHC-undiscovered pMSSM model sample can depend upon the prior used to generate 
the original model set. 



5 Note that the amount of fine-tuning in both our model sets was examined in some detail in our earlier 
work [7] 



36 





Figure 22: Two projections of the fine-tuning distributions for models in our flat (top panels) 
and log (bottom panels) prior sets. The top red histogram in all panels shows the result 
for the full model set while the subsequently lower histograms correspond to models not 
observed by the ATLAS E™ 1SS search analyses for various values of the integrated luminosity 
as indicated, assuming a background uncertainty of 50%. 
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4.2 Implications of pMSSM Searches for the Linear Collider 

If there are no clear SUSY signals as the 7 TeV LHC integrates more luminosity, the ques- 
tion arises whether the production (and study) of charged sparticles remains viable at the 
proposed 500 GeV Linear Collider (LC). Based on mSUGRA/CMSSM model coverage pro- 
jections from both the ATLAS [E] and CMS [29j Collaborations at 7 TeV with a 1 fb" 1 
integrated luminosity (as well as their results from the 2010 SUSY searches), it would seem 
very unlikely that either light sleptons or gauginos (other than perhaps the LSP itself) will 
remain kinematically accessible at a 500 GeV LC if nothing is observed. Here, we address 
the question whether this expectation also remains true for our pMSSM model sets. To 
this end we examine the set of flat and log prior models which are not detected by any of 
the ATLAS E^ 1SS search analyses for assumed values of both the integrated luminosity and 
background systematic error and then determine the part of the sparticle spectrum within 
these models which is kinematically accessible at a 500 GeV LCj^] A similar analysis can also 
be performed for a 1 TeV LC. We remind the reader that, e.g., in the flat prior case we have 
no sparticles heavier than ~ 1 TeV. 

We present the results of this analysis in various ways. Tables [7] and [8] show the 
number of sparticles of various species that are kinematically accessible at a 500 GeV or 
1 TeV LC within the subset of ATLAS-undetected models assuming C = 1(10) fb -1 with 
5B = 50(20)%, respectively. These two cases represent both a possibly conservative and 
a more optimistic performance for the LHC over the 2011-12 running period. Here we see 
several things: (i) In addition to the gauginos and sleptons, sparticles such as stops, sbottoms 
and other squarks are potentially almost as likely to also be kinematically accessible at the 
LC [32]. (ii) The number of models with kinematically accessible sparticles and their variety 
is significantly greater in the log prior sample as these models are more likely to have a 
lighter and more compressed sparticle spectrum (Hi) The difference between the two cases 
presented in these Tables is quite significant; in particular, we see that for the flat prior model 
set there is a huge depletion in the number of unobserved models at higher luminosity and 
with lower background systematics. There are extremely few flat prior models remaining at 
high luminosity with any accessible sparticles at a 500 GeV LC. (iv) Going from a 500 GeV 
to a 1 TeV LC substantially increases the number of models with kinematically accessible 
sparticles, especially in the flat prior case. It is clear that, at least for the flat prior model 
sample with luminosities in excess of 1 fb -1 at the LHC, that the 500 GeV LC does not seem 
to be a good place to study our pMSSM models if no signal for SUSY is found at the LHC 
in 2011-12. 

In order to study these LC results in more detail we examine their dependence on 
the LHC integrated luminosity; this is shown for the 500 GeV LC in Fig. [23] and for the 1 
TeV LC in Fig. [24| Here we display the fraction of the unobserved set of models that have 
a kinematically accessible sparticle of a particular variety. At a 500 GeV LC, this fraction 

6 Note that we have not performed any analysis here to ascertain whether or not a given kinematically 
accessible sparticle is actually observable at such a LC. See, however, the work in Refs. [3Q| f3Tj . 
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Table 7: Number of kinematically accessible sparticles from our set of 14623(1546) flat (log) 
prior pMSSM models that are unobservable by the ATLAS E™ ss searches assuming C = 1 
fb" 1 with SB = 50% for both a 500 GeV and 1 TeV LC. 
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Table 8: Same as the previous Table but now corresponding to the 672(663) undetected 
flat(log) prior models assuming C = 10 fb _1 with SB = 20%. 
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for charginos, stops, sbottoms, selectrons (or smuons) and staus in the flat prior model 
set is seen to decrease significantly as the LHC accumulates integrated luminosity without 
observation of a signal for SUSY. However, note that for the corresponding log prior model 
set, while the number of surviving models decreases with any corresponding increase in the 
LHC integrated luminosity (or with any decrease in the SM background uncertainty) as one 
would expect, the fraction of the surviving models with a kinematically accessible sparticle 
changes very little, if at all! For charginos at a 500 GeV LC this fraction is found to be quite 
large, ^ 75%, but it is found to be somewhat smaller for the other sparticles, e.g., ~ 25% 
for fi and ~ 7% for e L . 

One possible explanation of this unexpected behavior in the log prior sample is as 
follows: As we saw in the previous discussion of fine-tuning, in the log prior case, the removal 
of pMSSM models from the log prior set (as they are 'discovered' by ATLAS) must affect the 
various sparticle mass distributions in a roughly uniform manner. Otherwise the observed 
amount of fine-tuning would necessarily increase. However, in the flat prior case, models 
with lighter sparticles are preferentially 'discovered' by ATLAS searches. This hypothesis 
can qualitatively explain why there is no significant reduction in the fraction of the log prior 
models with kinematically accessible sparticles at a 500 GeV (or 1 TeV) LC (especially in the 
case of non-colored sparticles as we have seen above). It also simultaneously explains why 
the apparent amount of fine-tuning does not change appreciably as the LHC covers more of 
the log scan parameter space. Of course, as can be seen from these figures, at a 1 TeV LC 
a substantial portion of both the log and flat prior model sets which remain undiscovered 
at the LHC have sparticles which are kinematically accessible even at high LHC integrated 
luminosities. We further note that the fraction of LHC-unobserved models where no SUSY 
sparticles whatsoever are accessible at a LC (not even the LSP) is quite large for the flat 
model set at a 500 GeV LC, but is only in the 12-20% range for the corresponding log prior 
model set. On the other hand, at a 1 TeV LC, this value is seen to lie below ^ 1% for the 
log prior models, while for the flat prior set it remains below ~ 13%. Clearly, a 1 TeV LC 
will be far better at accessing the sparticles in our pMSSM model sets. 

In order to test this hypothesis, we show the mass distributions for the and \2 m 
both the flat and log prior model sets in Fig. 25|^ Indeed, we see that this hypothesis is true 



for non-colored states: in the flat prior set, models being observed by the ATLAS analyses 
mostly correspond to those with lower sparticle masses. On the other hand, for the log prior 
case, we see that the mass distributions for non-colored sparticles approximately maintain 
their overall shape as the models are observed, showing not much preference for the lighter 
sparticle masses. Fig. [26] shows, however, that in the case of colored sparticles, here for the 
gluinos and (one of the) squarks, which are most directly sensitive to most of the ATLAS 
E™ 1SS analyses, this same effect is somewhat less significant. In particular for the gluinos we 
see that even in the log prior case there is a significant loss in the fraction of models with 
lighter masses as the LHC integrated luminosity is increased. 



7 Note that there is nothing special about the choice of these two particular non-colored sparticles and 
the features that we will now describe are also found in the mass distributions of other sparticles. 
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Figure 23: Top left: Fractional number of undetected models with kinematically accessible 
xf at a 500 GeV LC as a function of the LHC integrated luminosity for flat (solid) and 
log(dashed) prior models. The red(green, blue) curves correspond to background systematic 
uncertainties of 50(100,20)%, respectively. Bottom left: Same as the previous panel but 
now for flat prior models only with 5B = 50% for the fi(blue), e# (green), e L (red), t\ (solid) 
and b\ (magenta). Bottom right: Same as the previous panel but now for log prior mod- 
els. Top right: Fraction of undetected models with no sparticles kinematically accessible. 
The magenta(green,blue) curves are for log prior models with SB = 20(50, 100)% while the 
blue (red, cyan) curves are the corresponding results for the flat prior set. 



42 




Luminosity (fb *) Luminosity (fb *) 



Figure 24: Same as in the previous Figure but now for a 1 TeV LC. 
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e L Mass Distribution for FLAT models failed for 50% error 



e L Mass Distribution for LOG models failed for 50% error 
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Figure 25: Mass distributions for eL(top) and ^(bottom) for the ATLAS-undetected 
flat(left) and log(right) prior models assuming 5B = 50% for different values of the LHC 
integrated luminosity as indicated. The top red histogram in each case corresponds to the 
original model sets before any of the ATLAS analyses are considered while the subsequently 
lower histograms correspond to those subsets of models undetected by the ATLAS E^ 1SS 
analyses at fixed integrated luminosities. 
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g Mass Distribution for FLAT models failed for 50% error 



g Mass Distribution for LOG models failed for 50% error 





u L Mass Distribution for FLAT models failed for 50% error 



u L Mass Distribution for LOG models failed for 50% error 
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Figure 26: Same as the previous Figure but now for gluinos and u^. 
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5 Summary and Conclusions 



In this paper we have analyzed the capability of the ATLAS £ , ™ lss -based SUSY analyses to 
discover Supersymmetry in a model independent fashion at the 7 TeV LHC. To this end, 
we tested these search channels on a large set of model points, ^ 71k, in the 19-dimensional 
parameter space of the pMSSM. This model sample contains a wide variety of properties and 
characteristics and provides a framework to explore the breadth of possible SUSY signatures 
at colliders and elsewhere. These models were generated in a previous work and comply with 
a set of minimal theoretical assumptions as well as the global precision electroweak, heavy 
flavor, collider searches, and astrophysical data sets. We simulated ten ATLAS E™ 1SS search 
channels, which were designed in the context of mSUGRA-based SUSY, and employed the 
SM backgrounds as provided directly by the ATLAS SUSY working group. We first checked 
that our analyses were in agreement with ATLAS results for the mSUGRA benchmark point 
that the collaboration had previously simulated. 

We passed our model sample through the ATLAS analysis chain and computed the 
significance of the signal for each model in each search channel. A significance of S > 5 
was used as the criteria for discovery in each channel; we employed the same numerical 
technique that ATLAS does for calculating this value. We found that the systematic error 
due to uncertainties in the size of the expected background made a substantial impact on 
model discovery for the range of expected integrated luminosities. In fact, some channels 
become systematics dominated at luminosities of order 5 — 10 fb _1 . Overall, for 1(10) fb _1 
of integrated luminosity roughly 80(95)% of the flat prior model sample is discoverable, 
assuming a 50% background systematic error. Larger (or smaller) systematic errors greatly 
reduce (or increase) this model coverage. We found that the 4j01 channel is the most powerful 
in terms of observing a signal, whereas the leptonic channels had a much reduced model 
coverage. This is due to the suppression of leptonic cascade decays appearing in our model 
sample compared to expectations from e.g., mSUGRA. The model coverage was worse for 
the log prior sample due to kinematic reasons. 

We explored the characteristics that caused a model to not be observed in these search 
channels. While production cross section values as related to the sparticle mass obviously 
plays a role, it does not tell the whole story. There are cases with low mass gluinos and/or 
squarks with large cross sections that are missed by these search analyses, while models 
with heavy masses and small cross sections are sometimes observed. We found that the 
mass splitting between the gluino/squarks and the LSP plays an important role in detecting 
models, and that this can sometimes be compensated by very large production rates or ISR. 
We also saw that subtleties in the sparticle spectrum can conspire to render a model to not 
be detected. A fraction of our model set contains detector-stable sparticles which appear at 
the end of their cascade decay chains and hence are not detected by the E^ lss -based searches. 
We studied the effectiveness of the planned stable charged particle searches in these cases 
and found that some, but not all, of these models will be discovered. 
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We briefly considered potential modifications to the ATLAS kinematic cuts in these 
E™ 1SS analyses that would improve their discovery potential. We studied the optimal cut on 
M e ff as well as for the of the leading jet and overall E^ 1SS . Our results indicate that the 
cuts for both the leading jet transverse momentum and the E^ 1SS could be increased from 
their nominal value without seriously impacting model coverage. 

Lastly, we studied the implications of a null result from the 7 TeV LHC run. We 
found that the degree of fine-tuning that would be imposed on the pMSSM depended on the 
choice of priors which generated the model sample, but overall would not be as large as in 
the case of mSUGRA. However, the expectations for sparticle production at a high energy 
Linear Collider would be greatly impacted if Supersymmetry is not discovered during this 
LHC run. Basic kinematics would essentially exclude sparticle production at a 500 GeV 
Linear Collider, and would point towards the need for a higher energy machine in order to 
study Supersymmetry. 

In summary, we find that the mSUGRA motivated £j? lss -based searches for Super- 
symmetry perform well over a larger and more complicated SUSY parameter space such as 
the pMSSM. However, there are some exceptions and coverage is not perfect. The details of 
the full sparticle spectrum play a very important role in the observability of a model. There 
are no blanket statements regarding the potential for discovery, or in setting a mass limit, 
that that can honestly be made. 

We anxiously await the discovery of Supersymmetry in the near future. 

Acknowledgments 

This project would not have been possible without the assistance and input from 
many people. 

The authors would like to thank the members of the ATLAS SUSY group, in par- 
ticular, S. Caron, P. de Jong, G. Polesello and G. Redlinger, for discussions and, most 
importantly, for providing us with additional information detailing the ATLAS-generated 
SM background distributions for their SUSY studies and the corresponding results for their 
mSUGRA benchmark analyses. 

We would also like to thank the SLAC ATLAS group for assistance with computing 
resources. 

We would like to thank T. Plehn for his help overcoming the special problems we had 
with the implementation of PROSPINO for our pMSSM model set. We would also like to 
thank A. Djoudai and J. Conway for similar help with the implementations of SDECAY and 
PGS, respectively. 

We would like to thank L. Dixon for discussions about the theoretical assessment of 
the size of the systematic errors for the SM backgrounds to SUSY signals. 

We would like to thank R. Cotta for important discussions related to our analysis. 



47 



We would like to thank J. Cogan for his assistance with setting up the ATLAS analysis 
code during the earlier part of this work. 

Work supported by the Department of Energy, Division of High Energy Physics, 
Contracts DE-AC02-76SF00515 DE-AC02-06CH11357, and DE-FG02-91ER40684, and by 
the BMBF "Verbundprojekt HEP-Theorie" under contract 05H09PDE. 



References 

[1] For recent overviews of LHC searches for new physics, see E. Gross (ATLAS) and H. 
Flaecher (CMS), talks given at LaThuile 2011, Les Rencont res de Physique de la Valle 
d'Aoste, 27 Feb-5 Mar, 2011. 



[2] V. Khachatryan et al [ CMS Collaboration ], |arXiv:1101.1628 [hep-ex]]; CMS Col- 



laboration, |arXiv: 1103.0953 [hep-ex]]; CMS Collaboration, | arXiv:1103.1348| [hep-ex]]; 
G. Aad et al. [ Atlas Collaboration ], |arXiv:1102.2357 [hep-ex]]; ATLAS Collaboration, 
[arXiv:1102.5290 [hep-ex]]. 



[3] H. E. Habe r and G. L. Kane, Phys. Rept. 117, 75 (1985); S. P. Martin, arXiv:hep-| 
ph/9709356; D. J. H. Chung, L. L. Everett, G. L. Kane, S. F. King, J. D. Lykken and 
L. T. Wang, Phys. Rept. 407, 1 (2005) |arX iv:hep-ph/0312378| ; L. Pape and D. Treille, 
Rept. Prog. Phys. 69, 2843 (2006); H. K. Dreiner, H. E. Haber, S. P. Martin, Phys. 
Rept. 494, 1-196 (2010); M. Drees, R. Godbole and P. Roy, Hackensack, USA: World 
Scientific (2004) 555 p; H. Baer and X. Tata, Cambridge, UK: Univ. Pr. (2006) 537 p. 

[4] G. Aad et al. [ The ATLAS Collaboration ], |arXiv:0901 .05121 [hep-ex]]. 

[5] G. L. Bayatian et al. [ CMS Collaboration ], J. Phys. G G34, 995-1579 (2007). 

[6] J. Alwall, M. -P. Le, M. Lisanti, J. G. Wacker, Phys. Lett. B666, 34-37 (2008). 
|arXiv:0803.0019 [hep-ph]]. 



[7] C. F. Berger, J. S. Gainer, J. L. Hewett et al, JHEP 0902, 023 (2009). |arXiv:0812. 0980 
[hep-ph]]. 



[8] J. Alwall, P. Schuster, N. Toro, Phys. Rev. D79, 075020 (2009). | arXiv: 0810.3921 [hep- 
ph]]; N. Arkani-Hamed, P. Schuster, N. Toro, J. Thaler, L. -T. Wang, B. Knuteson, 
S. Mrenna, [hep-ph/0703088 [HEP-PH]]. 

[9] J. A. Conley, J. S. Gainer, J. L. Hewett et al, [arXiv: 1009.2539 [hep-ph]]. 



[10] S. Akula, N. Chen, D. Feldman , M. Liu, Z. Liu, P. Nath, G. Peim, |arXiv: 1103. 1197 



[hep-ph]]; S. Heinemeyer, |arXiv: 1103.0 952 [hep-ph]]; B. C. Allanach, T. J. Khoo, 
C. G. Lester, S. L. Williams, [arXiv: 1103.0969] [hep-ph]]; P. Athron, S. F. King, 



48 



D. J. Miller, S. Moretti, R. Nevzorov, [arXiv: 1102.4363 [hep-ph]]; O. Buchmueller, 
R. Cavanaugh, D. Colling, A. De Roeck, M. J. Dolan, J. R. Ellis, H. Flacher, S. Heine- 
meyer et al, |arXiv:1102.4585 [hep-ph]]; S. Scopel, S. Choi, N. Fornengo, A. Bottino, 
|arXiv:1102.4033 [hep-ph]]; B. C. Allanach, [arXiv: 1102.3149 [hep-ph]]; D. Feldman, 
K. Freese, P. Nath, B. D. Nelson, G. Peim, | arXiv:1102.2548| [hep-ph]]; O. Buchmueller, 
R. Cavanaugh, D. Colling, A. De Roeck, M. J. Dolan, J. R. Ellis, H. Flacher, S. Heine- 
meyer et al, Eur. Phys. J. C71, 1583 (2011). |arXiv:1011.6Tl8| [hep-ph]]. 

[11] R. C. Cotta, J. A. Conley, J . S. Gainer, J. L. Hewett, T. G. Rizzo, JHEP 1101, 
064 (2011). | arXiv: 1007.55201 [hep-ph]]; R. C. Cotta, J. S. Gainer, J. L Hewett, 
T. G. Rizzo, Nucl. Phys. Proc. Suppl. 194, 133-144 (2009). [arXi v:0909.4088| [hep-ph]]; 
R. C. Cotta, J . S. Gainer, J. L. Hewett, T. G. Rizzo, New J. Phys. 11, 105026 (2009). 
|arXiv:0903.4"409] [hep-ph]] . 

[12] ATLAS Note, Prospects for Supersymmetry discovery based on inclusive searches at a 
7 TeV centre- of -mass energy with the ATLAS detector, Tech. Rep. ATL-PHYS-PUB- 
2010-010, CERN, Geneva (2010). 

[13] L. Pape, D. Treille, Rept. Prog. Phys. 69, 2843-3067 (2006); G. Weiglein et al. [ LHC/LC 
Study Group Collaboration ], Phys. Rept. 426, 47-358 (2006). [hep -ph/0410364| . 

[14] A. Djouadi, J. -L . Kneur, G. Moultaka, Comput. Phys. Commun. 176, 426-455 (2007). 
|hep-ph/0211331| . 

[15] G. DA mbrosio, G. F. G iudice, G. Isidori, A. Strumia, Nucl. Phys. B645, 155-187 
(2002). |hep-ph/020703 6|. 

[16] F. Abe et al. [ CDF Collaboration ] , Phys. Rev. D46, 1889-189 4 (1992); V . M. Abazov et 
al. [DO Collaboration], Phys. Rev. Lett. 102, 161802 (2009). |arXiv:0809]4472| [hep-ex]]. 

[17] V. M. Ab azov et al. [ DO Collaboration ], Phys. Lett. B660, 449-457 (2008). 
[arXiv:0712.3 805 [hep-ex]]. 

[18] We would like to thank the members of the ATLAS SUSY Physics Group for providing 
us with this information. 

[19] A. Djo uadi, M. M. Muhlleitner, M. Spira, Acta Phys. Polon. B38, 635-644 (2007). 
|hep-ph/0609292| . 

[20] W. Beenakker , R. Hopker, M. Spira, P. M. Zerwas, Nucl. Phys. B492, 51-103 (1997). 
|hep-ph/9610490[ ; W. Beenakker, M. Kramer, T. Plehn, M. Spira, P. M. Zerwas, Nucl. 
Phys. B515, 3-14 (1998). |hep-ph/971045 1|; W. Beenakker, M. Klasen, M. Kramer, 
T. Plehn, M. Spira, P. M. Zerwas, Phys! Rev. Lett. 83, 3780-3783 (1999). |liep- 
ph/9906298 [; M. Spira, |hep-ph/0211145| ; T. Plehn, Czech. J. Phys. 55, B213-B220 
(2005). |hep-ph/0410063|. 



49 



[21] P. M. Nadolsky, H. -L. Lai, Q. -H. Cao, J. Hust on, J. P umplin, D. Stump, W. -K. Tung, 
C. -P. Yuan, Phys. Rev. D78, 013004 (2008). |arXiv:0802.0007 [hep-ph]]. 



J. Conway, PGS4, http://www.physics.ucdavis.edu/~conway/research/software/pgs/pgs.html 



M. Fairbairn, A. C. Kraan, D. A. Milstead, T. Sjostrand, P. Z. Skands, T. Sloan, Phys. 
Rept. 438, 1-63 (2007). |hep-ph/0611040| . 

V. Khachatryan et al. [ CMS Collaboration ], | arXiv:1101.1645| [hep-ex]]. 

A. R. Raklev, Mod. Phys. Lett. A24, 1955-1969 (2009). [ arXiv:0908.0315| [hep-ph]]. 

I. Hinchliffe, F. E. Paige, M. D. Shapiro, J. Soderqvist and W. Yao, Phys. Rev. D 55, 
5520 (1997) [arXiv:hep-ph/9610544| . 



A. Strumia, | arXiv:1101.2195| [hep-ph]]. 

S. Cassel, D. M. Ghilencea, S. Kraml et al., |ar Xiv:1101.4664| [hep-ph]]. 
CMS Note 2010/008. 

C. F. Berger, J. S. Gainer, J. L. Hewett et al, Phys. Rev. D80, 095018 (2009). 
|arXiv:0712.2965] [hep-ph]] . 

C. F. Berger, J. S. Gainer, J. L. Hewett et a/., Phys. Lett. B677, 48-53 (2009). 
|arXiv:071 1.1374] [hep-ph]] . 

T. G. Rizzo, [arXiv: 1005.4626] [hep-ph]]. 



50 



